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Abstract. Recently there has been much interest in optomechanical devices for the 
■ production of macroscopic quantum states. Here we focus on a proposed scheme for 

. achieving macroscopic superpositions via nested interferometry. We consider the effects 

of finite temperature on the superposition produced. We also investigate in detail the 
scheme's feasibility for probing various novel decoherence mechanisms. 
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1. Introduction 

Optomechanical systems have long been investigated as a means of probing the quantum- 
to-classical transition in macroscopic devices [TJ [21 El HJ El El U\- However, it has 
generally proven difficult to meet all necessary conditions for such experiments. Firstly, 
a sideband-resolved device is required, allowing ground state cooling [SI El QUI Ell EE21 EE3] • 
Secondly, the device's coupling rate must be faster than the mechanical frequency [2J[B], 
in order to create a distinguishable state displaced by more than the device's zero point 
motion. Finally, the device must meet the strong coupling criterion, ensuring that 
single photons remain in the cavity long enough to cause significant effects [21 [TH [15] . 
In practice it is very difficult to meet all of these competing requirements simultaneously. 

The authors have recently proposed a method to create quantum superpositions 
in weakly coupled systems via postselected nested interferometry [16J. This method 
greatly relaxes the above requirements, allowing the creation of quantum superpositions 
with devices more easily in reach of current technology [13] , as well as possible tests of 
novel decoherence mechanisms 
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Here we consider the experimental requirements of the proposed nested 
interferometry scheme, investigating in detail its tolerance of finite temperature in the 
resonator and finite temperature in the surrounding environment. We also analyze 
in detail the time scale on which decoherence mechanisms operate, including both 
traditional environmentally induced decoherence [17] and proposed novel decoherence 
mechanisms [HI HH1 [201 ED E21 [23]. 



2. Nested interferometry 

Optomechanical systems evolve under the following Hamiltonian [2] 

"H = hhj a}a + hu m c^c — hga)a (c + c') , (1) 

with h defined as the reduced Planck's constant, oj the optical angular frequency, u m 
the mechanical angular frequency, the optomechanical coupling rate g = u q Xq/L, with 
the zero point motion xq = \Jh/ (2muu m ), a the optical annihilation operator, and c the 
mechanical annihilation operator. 

If a single photon is input to the cavity and the mechanical state begins in coherent 
state \'y) m , then the mechanical state will evolve as follows pQ: 

\m)rn =^ m \l(t)+a(t)) m (2) 

7(t) = ie ~ iuJmt 
a(t) = k(1 - Q~ lU)mt ) 
<p{t) = k 2 (u m t - sin u m t) , 

with k = g/oj m . Here we define the set of coherent states I7) as well as the single 
quantum-added coherent states (7, 1) [25] : 

00 n 

| 7 > -e-H 2 /^^|n} (3) 



Ov Tt 



ct| 7) exp(-|7P/2)^ ^ 

As detailed in [TB] the postselection is accomplished by means of an inner Mach- 
Zehnder interferometer (Fig. [T|). The single photon is input and split into both cavities 
by a beam splitter, creating state l/y/2(\l) a |0) fe + |0) a After weakly interacting 

(k ^ 1, a(t) 1) with the optomechanical resonator for time t, the state will be: 

\m) = j= |l) a |0> 6 \l(t) + a(t)) m + |0> a |1) 6 \j(t)) m } . (5) 



By postselecting for photons which exit the dark port, we select the l/y/2(\l) a |0) 



b 
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Figure 1. The inner interferometer is a Mach-Zehnder interferometer. The upper 
path contains Cavity A which has a weak optomechanical coupling to a resonator. 
In the absence of optomechanical interaction the interferometer is balanced and all 
light exits via the bright port. Postselecting only the photons which exit the normally 
dark port prepares the resonator in its excited state. Cavity B is used to match the 
spectrum and time delay of cavity A, and has no optomechanical interaction. 



|0) a |1) 6 ) component, and compute it to lowest order in k: 

= \ ¥ m lT(t) + a(t)) m - | 7 (f)> J 
« \ [e^fl(a(()) - l] h(t)) m 
« i [(1 - »k 7 sin( Wm t))(l + a(t)ct - a*(t)c) - l] | 7 (f)) m 

« 1 [« 7 (1 - cos^t)) | 7 (t)) m + a(t) VlTl 2 + 1 17(f), l) m ] (6) 

with D(rj) denned as the displacement operator. 

In the 7 = case, where the resonator has been cooled to its ground state, the 
above simplifies to a postselected state of: 

l«0) m = ^|l) m - (7) 

Thus in this case, the resonator is placed into the first excited state with probability 
|a;(£)| 2 /4. The weak interaction between the photon and the device is probabilistically 
amplified. 

2.1. Finite device temperature 

However, for a device of finite temperature, 7^0. Consider a mechanical resonator 
initially in a thermal state, a statistical mixture of coherent states: 

Pth = ^ / e-H 2 /^(| 7 >(7l)d 2 7 (8) 
7rn th J 
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where n t ^ is the average number of phonons: 

1 

and where fee represents the Boltzmann constant. Note that n t h is also the value of I7I 2 
averaged over the thermal distribution, Eqn. |HJ Note that in this subsection we will deal 
only with mechanical states and will thus drop the m subscript. 

For an initial coherent state, we will have created |"0 ps (*)) from Eqn. [6l a 
superposition between a small early component with mechanical state |^ ps (t)) and a 
large late component still in \"f(t)). 

To lowest order in k, the overall probability of successful postselection for an initial 
coherent state will be: 



(^ ps (t)|Vps(*)) f 

Note that | (-f(t)\ip ps (t)} 



[1 - cosu m t) 



1 2, 

+ 7T« 7 



2 sin 2 co m t 



(10) 



'l/4)/t 2 |7| 2 sin 2 ci; m t, precisely the second term of 
Eqn. [TOl Thus the first term represents our signal, while the second term represents a 
background noise of dark port events due to finite temperature rather than successfully 
conveying a phonon to the device. Averaging Eqn. [10] over the thermal distribution, 
Eqn. |H1 we arrive at: 

UJmt 1 



<<iMt)|iUf)>) 



'th 



2 • 2 

K sm 



-K 2 n t h sin 2 uj m t 
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So for the signal to be larger than the noise, we must have n t h <C 
4 [sin(u; m £/2)/ smu m t} 2 = sec 2 (u; m t/2). This implies that the nested interferometry pro- 
posal will only be successful if n <C 1, that is T <C hw m /kB- Thus, ground state cooling 
is essential for the success of this scheme. For a sideband-resolved device, this can be 
accomplished by driving the red (anti-Stokes) sideband of the cavity with a coherent 
beam M CB H2] • 



2.2. Nested interferometry 

The nested interferometry proposal [TB] aims to use this amplification to create 
macroscopic superposition states, doing so by means of the extended optical setup 
pictured in Fig. The postselection interferometer of Fig. [T] is nested in a larger 
interferometer with both an early and a late path. 

An experiment begins with the optomechanical device being cooled to its ground 
state by standard optomechanical cooling techniques [HI EE]- Single photons are input 
to the outer interferometer and are split into an early component and a late component. 
The late component enters the first delay line. The early component immediately enters 
the inner interferometer where it interacts with the device, and only the |1) component 
is passed through the dark port, entering a second equal length delay line. 

At this point, the late component is associated with mechanical component |0) m 
while the early component is associated with |l) m . These components are left to evolve 
freely for the length of the delay lines, which can, in principle, be arbitrarily long. During 
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Figure 2. The outer interferometer measures the coherence of superposition states by 
the use of matched time delays. The input pulse is split by polarizing beam splitters 
(PBS) into an early and late component each of which traverses the inner interferometer 
(see Fig. [1} . The early and late components are brought back together with a second 
delay line and the interference visibility is measured by varying the phase shift </>. 
During the interval between the early and late components the resonator will be in a 
(postsclected) superposition of excited and not-excited, and any decoherence during 
that time will reduce the final measured visibility. 



this time they may experience decoherence from either traditional environmentally 
induced decoherence [17] or one of many proposed novel decoherence mechanisms 

[igiia 123 El 1221 123|. 

Finally, the components exit the delay lines and the late component enters the 
interferometer. As before only the |1) component passes out of the dark port and we 
are left with both components in the |1) state, assuming no decoherence has occurred. 
At this point both components are interfered to check for visibility, allowing us to 
measure whether decoherence has taken place. 

This scheme has two advantages over previous schemes. First, it allows weakly 
coupled devices to be placed in superpositions by a single photon. Second, in principle, 
it allows observation of decoherence on an arbitrary time scale, as the delay lines can be 
varied. Previous schemes [21 El E] were limited in the time scales by both the mechanical 
period of oscillation and the cavity lifetime. This would require new devices to measure 
at different time scales. Though it may be difficult to determine the cause of the 
decoherence beyond any doubt, it will be possible to vary the temperature and the 
characteristics of the device, such as mass, frequency, mechanical quality factor, and 
optical finesse, allowing parameter dependence to be established. 



3. Decoherence 



Here we will review the various decoherence mechanisms to be considered in this paper. 
The devices to be considered are hypothetical optomechanical trampoline resonators 
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[T3"l IT6] , optimized for the nested interference scheme (Tab. [I]). 



Table 1. Wc include parameters for two trampoline resonators [13] close to being 
able to implement the scheme and two devices proposed in |16j that should allow 
the scheme to be implemented. The parameters are effective mass of the mechanical 
mode (ng). mechanical mode frequency (kHz), cavity length (cm), optical finesse of 
cavity, mechanical quality factor, environmentally induced decoherence temperature 
(K), k = g/u! m , and sideband-resolution measure uj m /T c . Proposed device no. 2 may 
be capable of observing novel decoherence mechanisms [2UJ UHl HH1 HE] ■ 



Device 


m 


fm 


L 


F 


Qrn 


Tetd 


K 




Tramp. #1 [T3] 


60 


158 


5 


38,000 


43,000 


0.3 


0.000034 


2.0 


Tramp. #2 [T3] 


110 


9.71 


5 


29,000 


940,000 


0.4 


0.0016 


0.09 


Proposed #1 [16] 


1 


300 


0.5 


300,000 


20,000 


0.3 


0.001 


3.0 


Proposed #2 p] 


100 


4.5 


5 


2,000,000 


2,000,000 


0.4 


0.005 


3.0 



3.1. Environmentally induced decoherence 

Most devices proposed for ground state cooling [HI [121 EE] require that the device be 
optically cooled below the temperature T env that the surrounding environment can reach 
by conventional cooling (there is one notable exception [27]). This is also true of the 
devices proposed in Tab. [1] 

In this situation, the mechanical resonator is modeled as coupled to an infinite bath 
of harmonic oscillators [TTJ, [8]. In the limit of ksT cm! ^> hu m , mechanical quality factor 
Q m ^> 1, and a Markovian regime with no memory effects in the bath, the bath degrees 
of freedom can be eliminated and the system can be described by the master equation 
for the reduced density matrix p [171 [28] : 



dt P h 



p,H T 



n [x,{p,p}]-^[x,[x,p}}, (12) 



h 1 ' L ^ ,rjJ h 2 

with H rcn the Hamiltonian from Eqn. [1] renormalized by the interaction of the device 
and the bath, the damping coefficient 7 m = oo m /Q m , and the diffusion coefficient 
D = 2m~f m k-BT em . The first term represents the unitary evolution of the system 
under the Hamiltonian from Eqn. [H while the second term represents the damping 
and the third term represents the diffusion. In the macroscopic regime the diffusion 
term proportional to D/h 2 dominates Eqn. [12] [TTJ [8]. Thus the resulting time scale for 
decoherence is: 

TEID ~ D{Kxf ~ 2^' (13) 

with the superposition size Ax = Xq. It is helpful at this point to define an 

environmentally induced decoherence temperature [S]: 

hu rn Q m , . 

Jeid — — ; • 
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We note that the inverse of the decoherence time scale is r^ D = 2uj m {T cxw /T-Em). Thus 
for the environmentally induced decoherence to act on a time scale slower than the 
mechanical frequency it is necessary that T cnv <C Ieid- 

We will consider EID with a base temperature of T env = 1 mK, obtainable with 
a dilution refrigerator. For this case, for the 300 kHz device, teid ~ 150 /is. For the 
4.5 kHz device, teid ~ 15 ms. 



3.2. Gravitationally induced decoherence 

Gravitationally induced decoherence, proposed independently by Diosi [18] and Penrose 
[T9] . is a type of decoherence caused by an object in superposition's perturbation of 
spacetime. The time scale for such decoherence is: 

r P = h/Ap (15) 

with the Ap defined as follows: 

A P = 47TC [f (Pl(^)-P 2 (y(Pl(^)-P2(f)) d 3 xd 3 y; (16) 

J J \x-y\ 
with pi(x) and p 2 (x) the mass distributions of the two superposed states. 

As in [8], we model the system as set of spheres representing nuclei. The 
Penrose energy for one sphere is given by A p = 47r(.E' 12 + — E® l — E^), with 
Emn = ~G II Pm(x)p n (x)/\x — y\d 3 xd 3 y. The spheres considered are far enough apart 
and displaced little enough that their most significant interaction is with themselves, 
and not neighboring spheres. This means that we can merely multiply by the number of 
spheres, M/m, to get the total energy A P = (M/m)A P = An(Ei y2 + E 2j i — E^i — £'2,2), 
with £ m , n = (M/m)E^ n 

For all cases, we will consider two spherical mass distributions with radii a 
equal to the size of the specific mass distribution that will be chosen, separated by 
Ax = Xq = \/h/ (2muj), the zero point motion of the resonator. Note that this is 
mathematically equivalent to the model of one sphere at x = for |0) m , and two half- 
mass spheres at x = ±xo for \l) m - 

As the radius of the two spheres will be greater than xq regardless of mass 
distribution used, there will always be significant overlap in the distributions. This 
will greatly complicate evaluation of Eqn. [T6j This has no effect on the self-energy 
terms but does affect the interaction terms. The 1/r potential between overlapping 
spheres has been evaluated previously [29J: 

GMm/Ax if Ax > 2a, 

12a 2 - 5Ax 2 Ax 5 - 30Ax 3 a 2 



-GMm 



(17) 

if < Ax < 2a. 



10a 3 160a 6 

For the E^i and E 2y 2 terms, we can just plug Ax = into Eqn. [T71 This gives 

El,l ~ ^2,2 —■ 

There is considerable theoretical disagreement about the proper mass distribution 
to use for gravitationally induced decoherence J2J EH1 El EH E2] • Previous papers have 



Macroscopic superpositions via nested interferometry 



8 



used the zero point motion of the resonator itself, the nuclear radius of the nuclei making 
up the resonator, the zero point motion of the nuclei making up the resonator, and a 
completely homogeneous mass with no nuclear granularity. At this point, we will define 
the mass distributions to be considered in this paper. 

3.2.1. Zero point motion of resonator Zero point motion is defined as: 

a = Xo = \l^£ (18) 
For the 300 kHz device, a = 5.3 fm. For the 4.5 kHz device, a = 4.3 fm. 

For this case, for the 300 kHz device, rp ~ 3.5 ms. For the 4.5 kHz device, rp ~ 28 
/is. This type of decoherence might potentially be testable in the 4.5 kHz device, as it 
is faster than EID. 



3. 2. 2. Radius of tantalum The atomic nucleus has a size of approximately [33] : 

a = r a A l l\ (19) 

with r = 1.25 fm and A the atomic mass number. Since the largest component of the 
mass of a Ta20s/Si02 dielectric mirror will be tantalum, we will make the simplifying 
assumption that the mirrors are composed of tantalum. For tantalum, A = 181, so 
a 7 fm. 

For this case, for the 300 kHz device, rp ps 7.1 ms. For the 4.5 kHz device, rp ~ 100 
/is. This type of decoherence might potentially be testable in the 4.5 kHz device, as it 
is faster than EID. 



3.2.3. Zero point motion of nuclei In the Debye model, the zero point motion of nuclei 
in a lattice is given (Eqn. 12.3.10 in [34]): 

3h 

a = Xo >- = «¥• (20) 
with G d the Debye temperature and M the atomic mass. Since the largest component of 
the mass of a Ta 2 05/Si02 dielectric mirror will be tantalum, we will make the simplifying 
assumption that the mirrors are composed of tantalum. The Debye temperature of 
tantalum is G^ = 240 K [35], and the atomic mass M = 181 amu. Thus a « 5 pm. 

For this case, for the 300 kHz device, rp ps 1.8 x 10 6 s. For the 4.5 kHz device, 
rp w 28 x 10 3 s. This type of decoherence would not be testable, as it is slower than 
EID in both devices. 



3.2.4- Homogeneous mass Some have even proposed modeling the resonator as a 
perfectly homogeneous mass with no nuclear granularity [311 [32]. In general this sets an 
extremely high bar for the decoherence times, but we will compute it for completeness. 
In this case we will model the mass as a single sphere of radius a = 30 /im (compared 
to a 60 /iin diameter cylinder) with mass 60 ng. It is as though the mirror is composed 
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of one very large nucleus. Though the shape is not correct, this model will suffice for 
an order of magnitude estimate. This can be represented by setting the nuclear mass m 
equal to the resonator mass M in Eqn. [T71 

For this case, for the 300 kHz device, rp pa 12 x 10 9 s. For the 4.5 kHz device, 
Tp ~ 1.8 x 10 12 s. This type of decoherence would not be testable, as it is slower than 
EID in both devices. 



3.3. Continuous Spontaneous Localization 

Continuous spontaneous localization is a proposed position-localized decoherence 
mechanism in which a nonlinear stochastic classical field interacts with objects causing 
collapse of macroscopic superpositions. Proposed by Ghirardi, Rimini, Weber and Pearle 
[221 123] . the master equation and decay rate for position-localized decoherence have the 
following form [321 HH [221 [231 HH EU [261 12Q] : 

d % 

(a; | p \x') = — (x\ [p, H] \x') — T(x — x') (x\ p \x') (21) 



dt N 1 r 1 ' h 
r(x) =7 



exp 



x 2 



(22) 
(23) 



4a 2 , 

Ax 2 if x < 2a, 
7 if x 3> 2a. 

with T(x) the decay rate, A = 7/(4a 2 ) the localization parameter, 7 the localization 
strength, and a the localization distance. In all cases, the trampoline resonators 
considered are in the x -C 2a limit. For the single nucleon case, the continuous 
spontaneous localization model |23] gives values ocsl = 100 nm and 7^sl = 

10- 16 Hz 

based on phenomenological arguments. 

Following [361 E2], the value of the localization parameter A C sl can be shown to 

be: 

Acsl—t^/^M (24) 

with M the resonator mass, mo the nucleon mass, R the radius of the sphere and 
f(R,b,a) a parameter depending on the geometry of the device. Disk geometry was 
considered in [36]. For motion perpendicular to the disk face / is evaluated (see [36] . 
Sec. 5.2, App. A, and Eqn. A. 11): 

/o„\4/ 9 \2 r R/2a r R/2a 

f(R,b,a)=4(^j fyj [1 - e- 62 / 4 " 2 ] J xdx J x'dx'e^^ I (2xx') (25) 

with R the disk radius, b the disk thickness, Iq{x) the n = modified Bessel function 
of the first kind, and a the localization distance (for CSL, acsL = 100 nm). In the 
(R/2a) 2 > 1 and (6/2a) 2 > 1 limits, applicable in this case, / pa (2a/R) 2 (2a/b) 2 . 

Thus, for the 300 kHz device, using a thickness of ~ 5 /im and a radius of ~ 4 urn 
(values consistent with the proposed finesse and mass), we obtain a decoherence time 
of order tcsl = 10 7 s. For the 4.5 kHz device, using a thickness of ~ 5 pm and a radius 
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of ~ 40 /an, we obtain a decoherence time of order tcsl = 1.5 x 10 5 s. This type of 
decoherence would not be testable, as it is slower than EID in both devices. 



3.4- Quantum gravity 

It has been proposed that quantum gravity might cause a form of position-localized 
decoherence due to coupling of the system to spacetime foam. This was first proposed 
by Ellis, Nanopoulos, Hagelin, and Srednicki [20] and subsequently elaborated j2TJ [26] 
with others. Notably, this model is phenomenologically equivalent to the CSL model 
with altered values for the constants (32]: oqq = hm-p / 2cm^ with mp = ^hc/G the 
Planck mass, and 7q G = Aa^ G c i mQ/ti?w?p. This gives us: 

A QG = 2 . 2 f{R, b, a) = f(R, b, a) (26) 



ml 4a 2 QG h 3 



m 



with f(R, b, a) as in Eqn. [251 However, since R <C Oqg and b <C Oqg, we can set / to 1 

AQG ~ ~hhnY (27) 
Thus, for the 300 kHz device, using a thickness of ~ 5 /im and a radius of ~ 4 /im, 
we get a decoherence time of order tqg = 7.1 s. For the 4.5 kHz device, using a thickness 
of ~ 5 /im and a radius of ~ 40 /im, we get a decoherence time of order tqg = 1.1 ms. 
This type of decoherence might potentially be testable in the 4.5 kHz device, as it is 
faster than EID. 



4. Conclusion 



In conclusion, we have presented an analysis of the experimental requirements of 
the nested interferometry scheme [16]. The scheme allows for the creation of 
macroscopic superpositions in weakly coupled systems, and allows for investigation of 
their decoherence on arbitrary time scales limited only by external delay lines. In 
particular, we investigate the temperature dependence of the scheme and find that 
ground state cooling is necessary for implementation. We also investigate the time 
scales on which proposed novel decoherence mechanisms would be expected to operate. 
We conclude that two proposed versions of gravitationally induced decoherence [HI [19] 
are testable, and that quantum gravitational decoherence [201 EI] is testable by this 
scheme. 
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